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In practice, physical membranes are exposed to a certain amount of external strain (tension 
or compression), due to the environment where they are placed. As a result, the behavior of 
the phonon modes of the membrane is modified. We show that anharmonic effects in stiff two- 
dimensional membranes are highly suppressed under the application of tension. For this, we consider 
the anharmonic coupling between bending and stretching modes in the self-consistent screening 
approximation (SCSA), and compare the obtained height-height correlation function in the SCSA 
to the corresponding harmonic propagator. The elasticity theory results are compared to atomistic 
Monte Carlo simulations for a graphene membrane under tension. We find that, while rather high 
values of strain are needed to avoid anharmonicity in soft membranes, strain fields less than 1% are 
enough to suppress all the anharmonic effects in stiff membranes, as graphene. 

PACS numbers: 81.05.ue, 68.60.Dv, 63.20.Ry, 46.70.Hg 



I. INTRODUCTION 

The study of mechanical properties of physical mem- 
branes, which are two-dimensional surfaces embedded in 
three-dimensional space, has a prominent experimental 
platform in graphenepJ a single layer of c arb on atoms ar- 
ranged in a hexagonal crystalline orderP^l Graphene is 
a stiff membrane, which presents long wavelength modu- 
lations of the out-of-plane displacements, commonly re- 
ferred as ripplesP® The impact of corrugation on elec- 
tronic transport, as well as the mechanical properties of 
graphene, are subjects of intense investigation!-^^ 

Of special interest is to understand the effects of an 
external strain applied to the membrane.^ This is so be- 
cause most of the graphene samples are subject to some 
finite amount of strain, due either to the pinning to the 
substrate (for samples on Si02, for example) or to the 
electrostatic force due to the gate on suspended sam- 
ples. In particular, the existe nce o f tension affects the 
dispersion of flexural phonons! ' ' In fact, whereas in 
the harmonic approximation and in the absence of strain, 
the dispersion relation of flexural phonons is quadratic, 
ujfi(q) ~ q 2 , strain introduces a characteristic wave- 
vector q s where the dispersion changes from linear (for 
q < q s ) to quadratic (for q > q s ). However, anharmonic 
coupling between bending and stretching modes is impor- 
tant and leads to a further renormalization of the mode 
dispersion, especially at long wavelengths .EHUD 

In this paper, we study the effect of tension on the 
flexural phonons of a 2D membrane. For this aim, we 
include a strain field in the free energy, which is stud- 
ied first in the harmonic approximation and then includ- 
ing anharmonic eff ects, u sing the self-consistent screen- 
ing approximation^^ 21 (SCSA). The results of a stiff 
membrane as graphene are compared to those for a softer 
membrane, for which the bending rigidity has been highly 
reduced. The validity of the continuum elastic the- 
ory is checked by comparing the SCSA results to atom- 
istic Monte Carlo (MC) simulations. Our numerical re- 



sults show that, for stiff membranes as graphene, small 
amounts of tension can be used to suppress the anhar- 
monic effects. The case of compressional strain is much 
more complex due to its highly non-linear behavior,EHHl 
and it cannot be accounted for in the SCSA. However, 
we can still use atomistic MC simulations for this case 
and, in fact, we find a highly non-trivial behavior for 
the correlation function of a compressed graphene mem- 
brane, with no crossover to a power-law behavior as in 
the tensioned case. 

The paper is organized as follows. In Sec. |H]we dis- 
cuss, in the harmonic approximation, the effect of an 
external strain in the system, and compare the results to 
the unstrained case. In Sec. Mil we consider the anhar- 
monic coupling between bending and stretching modes 
in the SCSA. In Sec. IIVI the results of the continuum 
elastic theory are compared to atomistic MC simulations 
for the height-height correlation function. The MC sim- 
ulation for a compressed membrane is also included. The 
main conclusions of our work are summarized in Sec. IVl 



II. HARMONIC APPROXIMATION 

In the absence of any external strain, the flat phase 
of a 2D membrane at sufficiently long scales is well de- 
scribed by a free energy that is a sum of a bending and 
a stretching parlP3 

J"[u, h] = \J d 2 v [k (V 2 h) 2 + 2^u% + Xu 2 aa ] (1) 

where k is the bending rigidity, A and [i are the first Lame 
constant and the shear modulus, respectively,^ and u a p 
is the internal strain tensor 

u a p ~ ^{d a Ufs + dpu a + d a hdph). (2) 

In the harmonic approximation, the bending and stretch- 
ing modes are decoupled. From Eq. (JlJ, one can calcu- 
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late the correlation function for the out-of-plane displace- 
ments h(r) which, in Fourier space, reads 



(IMq)l 2 )«=o 



k B T 



(3) 



where k B is the Boltzman constant, T is the temperature, 
and the suffix u — in the average denotes the absence 
of any external strain. The effect of an external strain 
applied to the membrane is modeled by the inclusion of a 
new term in the theory, r a p, that couples to the internal 
strain tensor. Therefore, we use the following expansion 
for the free energy 
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where r Q/ 3 = \5 a pu e *p + 2/xu^' is expressed in terms 
of the external strain tensor u^g. In the harmonic ap- 
proximation, the Fourier component of the height-height 
correlation function reads simply 



G (q) = (\h( q )\ 2 ) u 



q 2 («g a + A«»* + 2^^ 



(5) 

For an isotropic expansion and in the long wavelength 
limit, we can approximate 



l af3 



= u5 a f 
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where u = SS/2S accounts for the uniform dilation of the 
membrane, where S is the membrane surface and SS is 
the change in area due to the application of strain. This 
reduces Go(q) to 



Go(q) = 



k B T 



q 2 [nq 2 + 2(A + 



(7) 



where r = 2(A + fi)u is the stress of the system. No- 
tice that for the unstrained case (u — 0), as given by 
Eq. ([3]), the mean square amplitude of the out-of-plane 
displacement diverges, in the harmonic approximation, 
asP3 (h 2 ) u=0 oc L 2 , where (h 2 ) = £ q (IMq)| 2 ) and L is 
the sample size. Furthermore, the normal-normal cor- 
relation (n(r) • n(0)) diverges logarithmically as r — > 00. 
However, for a membrane under uniform dilation, we can 
use Eq. ([7]) and obtain 
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where q m in — 2tt / L is an infrared cutoff of the order of 
the inverse sample size L. Furthermore, from the normal- 
normal correlation function (|n(q)| 2 ) = kBT/[n(q 2 +q 2 )], 
where q u = \/2(A + Jtju/K we obtain that 

k R T f d 2 q e iqr k B T 



(n(r) • n(0)) 



K J (2tt) 2 q 2 + ql 
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FIG. 1. (Color online) Correlation function in the harmonic 
approximation, Go(q) (blue dashed line) and renormalized 
correlation function in the SCSA G(q) (full red line) for sev- 
eral values of strain u from to 10~ 2 . Dotted-dashed vertical 
lines indicate q c ~ 0.24 A -1 , according to the Ginzburg cri- 
terion, Eq. l |16[ ). Dotted vertical lines indicates the position 
of g^ 1 and dashed vertical lines point the position of q B (see 
text). 



where Kq(x) is a modified Bessel function of the second 
kind. Taking into account that, for x ^> 1, Kq(x) ~ 
\J n/2xe~ x , we obtain that, for r — > 00 



(n(r) • n(0)) 



k R T 
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K \/8nq u r 
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Therefore, the application of external strain, as expected, 
guarantees the long range 2D order of the membrane. 
Furthermore, q~ 1 defines a length-scale that separates 
the strain dominated from the unstrained regions of the 
correlation functions. 



III. ANHARMONIC EFFECTS: CONTINUUM 
ELASTIC THEORY IN THE SCSA 

In the previous section we have seen that the applica- 
tion of an external strain stabilizes, even in the harmonic 
approximation, the flat phase of a 2D membrane. But 
even in the absence of strain, it is known that the flat 
phase is stable. This is due to the anharmonic coupling 
between bending and stretching modesP-2 Therefore, an- 
harmonic effects lead to a further rcnormalization of the 
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characteristic lengths and elastic constants discussed in 
the previous section. In the following, we study the ef- 
fect of anharmonicity in the correlation function of a 
strained 2D membrane. First, one notices that the in- 
plane phonons in the free energy Eq. Q can be inte- 
grated out exactly, what allows us to write an effective 



action in terms only of the h fields^ 

1 



Feff[h,T a p] = J * d 2 



K (V 2 ft.) + Tafidahdph 



l -Y J d 2 r(P^d a hd h + P^u^) 2 



(11) 

where Y = 4/i(/x + A)/ (2fi + A) is the 2D Young modulus 
and Pj^ = 5 a p—d a dp/S/ 2 is the transverse projection op- 
erator. Eq. ( 11 ) can be expressed in terms of the Fourier 



components of the height field, h(q). To the lowest order 
in u ext we obtain 



Feff[h,T aP ] = \ J ^k 2 [nk 2 + \u% + 2^</) \h(V\ 2 

+ l Y Jw^JjSf 2 'J ^^^(q) fc i«^fc3 7 ^[^(ki)/ i (k 2 )][/ i (k 3 )/ i (k 4 )] 



(12) 



where q = ki + k 2 and in the second term ki + k 2 + 
k 3 + k 4 = 0. The first line of Eq. ( 12 1 is nothing but the 



bending part of the free energy in the harmonic approxi- 
mation, from which we have defined the non-interacting 
correlation function Go(q), as in Eq. ^ of the previous 
section. The second and third lines build the interac- 
tion term of the theory. The second line of Eq. ( 12 ) 



accounts for the four point vertex, whereas the last term 
of this equation leads to a two-point vertex that renor- 
malize the propagator. This problem is similar to that of 
a polymerized membrane with long-range disorder, and 
can be treated in the SCSAP 

For stiff membranes as graphene, the anharmonic ef- 
fects are quickly suppressed under the application of 
strain. Therefore, we assume that the renormalization of 
the propagator due to the vertices associated to the last 
term of Eq. ( 12 ) is weak, and that we can neglect this 



class of diagrams in our calculations. We will see in Sec. 
IIVI that the correlation functions calculated with this as- 
sumption in the SCSA agree well with those obtained 
from atomistic MC simulations, justifying the simplifica- 
tion. Then, the renormalized correlation function can be 
calculated from a closed self-consiste nt set of two coupled 
integral equations for the self-energjff^HH 

f d 2 q ~ 

S(k) = 2k a k [j k 1 k & J j^-^R a0 ^ s (q)G(k.- q) (13) 

(14) 

where G _1 (q) = G^q) + S(q) is the inverse of the 
dressed propagator, H a j3 tl s(q) are the vacuum polariza- 



tion functions, 



rf 2 k 
(27T) 2 



k a k /3 (k 7 -q 7 )(k s -qs)G(k)G(q-'k), 

(15) 
- P T P T ) 



= (Y/2-ripT p pT+( P /2)(P^Pj 5 + P'+ s P'P 
is the unrenormalized four-point interaction vertex and 
Rap, 75 (q) is the screened interaction. 

The set of equations ( 13 1-( 14 ) can be solved, at any 



wave- vector, following the method introduced in Ref. [221 
Fig. [I] shows the momentum dependence of Go(q) and 
G(q) for different values of strain, u = 0, 10~ 2 . First, 
one notices in Fig. [IJa) that, without strain, the har- 
monic approximation is valid only in the short wave- 
length region, where G(q) ps Go (q).^ The Ginzburg cri- 
terion, which considers only the first order correction to 
the correlation function, allows to estimate the character- 
istic wave-vector q c above which the harmonic behavior 
applies. For 2D membranes, q c is approximately given 



3k B TY 
8ttk 2 ' 



(16) 



For the parameters of graphene, q c ~ 0.24 A -1 at room 
temperature, a value which is shown by the vertical dot- 
dashed lines in Fig. [TJ This is the characteristic wave- 
vector at which the renormalized correlation function 
G(q) (full red line) separates from the harmonic approx- 
imation Go(q) (dashed blue line), pointing out that an- 
harmonic effects are very important at long length scales. 
If we consider the effect of external strain on the mem- 
brane, we still can distinguish between the harmonic and 
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The temperature dependence of the parameter A of 
obtained from the SCSA correlation functions at 
different temperatures (red dots). The dashed line is a fitting 
to Eq. (22 1. For these plots, we have used the wave- vector 
q = 1CT 



A 



the anharmonic regimes, as shown in Fig. [T|b)-(f)!^ 
However, there also exists a characteristic scale at which 
the behavior of flexural phonons is dominated by strain 
effects. This scale manifests itself as a change in the slope 
of the correlation functions: Go(q) changes from ~ q~ 4 
to q~ 2 and G(q) changes from q~ 4+T > also to q~ 2 , where 
r\ is a characteristic exponent. In fact, we can observe 
from Fig. [ljb)-(f) how the region of intermediate mo- 
menta where the two functions Go(q) and G(q) are dif- 
ferent, and therefore anharmonic effects are important, 
is reduced as the value of u grows. From these results 
we see that very small amounts of strain are enough to 
suppress the anharmonic effects. 

These results can be used to study the effect of strain 
on flexural (out-of-plane) phonons. The dispersion re- 
lation for flexural phonons of a 2D membrane under 
isotropic tension can be written as 



Mq) 2(A + n) 

W /Kq) = y <r + u — - — i 



(17) 



where p is the density and «(q) is the bending rigidity. 
In the harmonic approximation, ft(q) = n and the dis- 
persion changes from linear to quadratic at a wave- vector 
equal to 



2u(fi + A) 



(18) 



This characteristic wave-vector is denoted by the verti- 
cal dotted lines in Fig. [ljb)-(f)p2l However, anharmonic 
effects are important at long scales. To obtain analyti- 
cal results, we use the effective Dyson equation for the 
correlation functiorP22l 



(19) 



G- 1 (q) = G - 1 (q) + S a (q), 



where G a (q) is an approximated correlation function 
dressed by the self-energy E a (q), which is approximated 



by 



S„(q) = V 



(20) 



where A is some numerical factor, r/ « 0.82p^j and qo = 
2TTy/Y/n. From the approximation Eq. (19) one can 
obtain the renormalized bending rigidity 



Kfl(q) 



k B TA 



(21) 



It is important to mention that the coefficient A is tem- 
perature dependent. Notice that anharmonic effects are 
present in G Q (q) below a characteristic wave- vector q* , 
which is solution of T, a (q*) « Gq 1 (</*). Assuming that q c 
is the only crossover wave- vector from harmonic to anhar- 
monic behavior, and that q* ~ q c , then one can easily ob- 
tain th at th e temperature dependence of the parameter A 
in Eq. ( 20 ) follows the power-law A cx (fc B T/K)? -1 PUBy 
fitting the SCSA correlation function for different tem- 
peratures to Eq. (19), we find the dependence of the 



parameter A on temperature, and the results are shown 
in Fig. [2] (red dots) . This allows to define an approxi- 
mate expression for the adimensional parameter A, which 
is (using the elastic constants valid for graphene) 



A 4.6T[KY- 



(22) 



where T[K] is the temperature expressed in Kelvin. The 
results are shown in Fig. [2] by the dashed line, which fits 
rather well the values of A obtained numerically. This 
confirms that the assumption q* ~ q c is indeed valid 
within the SCSA. The main message is that the bending 
rigidity Eq. (21) grows with temperature as 



T r,/2_ 



(23) 



This power-law behavior is similar to the temperature 
dependence f ound by Monte Carlo simulations in the har- 
monic regime PSH However we emphasize that here we are 
assuming that the parameters k, /i and A of the Hamilto- 
nian ( lj are independent of temperature. However, while 
A and /x are only weakly dependent on T, the tempera- 
ture dependence of the bending rigidity n found in MC 
simulations is rather strongpSl and this is not accounted 
for by Eq. ( p3| . The origin of this T-dependence is prob- 
ably beyond the continuum medium approximation, and 
it lies beyond the scope of this work. 



Then, from Eq. (|2l|) one observes that the slope of the 

2-77/2 



dispersion relation cjj/(q) changes from ~ q to 
at the wave-vector solution of 
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Qo 



q 2 s =2u(X + ^). 



(24) 



The values of q s for the values of strains studied here 
are shown by the vertical dashed lines in Fig. [ljb)-(f ) . 
Notice that the characteristic wave-vector obtained by 
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the approximation to the bending rigidity Eq. (21 1 agrees 
well with the exact result of the SCSA equations. 

From the previous expressions, and imposing q s = q c , 
it is possible to find the critical value for the strain that 
is enough to suppress the anharmonic effects completely, 
at any wave-vector, and this is 



3k B T 



47T k(2^ + A)' 



(25) 



For the parameters of graphene at room temperature, 
this corresponds to u c s» 0.0025. In fact, notice that q^ 
and q s already coincide for u = 10~ 2 > u c and that both 
are to the right of q c (q s ,q% > q c ) [Fig. [ljf )] , pointing 
out that anharmonic effects are already absent for this 
value (~ 1%) of external strain. 

Finally, we compare the results for graphene to those 
of a softer membrane. In Fig. [3]Ja) we show Go(q) and 
G(q) for graphene, and for a membrane with the same fi 
and A as graphene, but with a bending rigidity k which 
is 1% of the corresponding for graphene. In this case, 
we see that q s ° ft > q9 ra P h , as seen by the position of 
the vertical dotted-dashed lines (red for graphene and 
black for the soft membrane). This means that anhar- 
monic effects manifest themselves at larger wave-vectors 
for a soft membrane. Furthermore, the change in slope 
of the harmonic Go(q) occurs at higher wave- vectors for 
the soft membrane (vertical black dotted line) as com- 
pared to graphene (vertical red dotted line). However, 
the out-of-plane component of the dispersion for flexural 
phonons dominates in a wider region of momenta for the 



soft membrane as compared to graphene, q 



soft 



graph 



as it can be seen by the relative position of q s for graphene 
(dashed red line) with respect to that of a soft membrane 
(dashed black line). Notice that in the latter case, the 
strain necessary to suppress all the anharmonic effects is 
u c ss 0.25, also two orders of magnitude larger than for 
graphene. In Fig. |3jb) we compare the correlation func- 
tions of graphene to those of a softer membrane where 
not only the bending rigidity k, but also the Lame con- 
stants A and /i have been reduced to 1% of their values in 
graphene. The situation is similar to that described for 
Fig. [3|a), with the difference that the wave- vectors at 
which anharmonic and strain effects appear are reduced, 
as it can be seen by the respective shifts to the left of the 
black dotted-dashed and dashed lines in Fig. [3](b) with 
respect to (a). Furthermore, for the parameters of Fig. 
^b), the change in slope of the harmonic Go(q) is the 
same in the two cases (as shown by the vertical dotted 
line). Finally, we notice that increasing the temperature 
acts like an effective softening of the membrane, due to 
the reduction of the ratio n/kgT. 



IV. COMPARISON TO ATOMISTIC MONTE 
CARLO SIMULATIONS 

In this section we compare the results obtained by us- 
ing the continuum elastic theory methods as described 
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FIG. 3. (Color online) Go(q) for graphene (dashed blue line) 
and for a softer membrane (dashed gray line), and G(q) for 
graphene (full red line) and for a softer membrane (full black 
line). In the two cases, u = 10 -8 . In (a) we have used, for the 
soft membrane, 1/100 times the bending rigidity k valid for 
graphene at this temperature, whereas {i and A are the same 
as in graphene. In this case, q c ~ 0.24 A -1 for graphene 
(vertical red dotted-dashed line) and q c ~ 24 A -1 for the 
soft membrane (vertical black dotted-dashed line). In (b) we 
compare the correlation functions for graphene to a softer 
membrane for which all the elastic constants are reduced to 
a 1% their value in graphene. In this case, q c m 2.4 A -1 
for the soft membrane, as indicated by the position of the 
vertical black dotted-dashed line. The vertical dotted and 
dashed lines represent the positions of and q s respectively, 
as in Fig. [l] 



in Sec. [TT| and |III| with the results of the Monte Carlo 
simulations of graphene. The correlation function G(q) 
for graphene has been calculated as described in Ref. 
1191 for unstrained graphene by means of MC simulations 
based on an accurate interatomic potential for carbonP^ 
The simulations are done for a sample of 37888 atoms 
in a roughly square sample of 314.82 x 315.24 A 2 , in 
the NPT isothermal isobaric ensemble. The simula- 
tions for the strained case were done for smaller sam- 
ples of 8640 atoms (147.57 x 153.36 A 2 ), which limits 
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FIG. 4. (Color online) Comparison of the normal-normal cor- 
relation function (jn(q)| 2 ) = g 2 G(q) obtained f rom contin- 
uum elastic theory, as described in Sec. [EI] and [Til] (dashed 
lines), to atomistic MC simulations (solid lines), for different 
values of external strain. 



FIG. 5. (Color online) Comparison of the normal-normal cor- 
relation function (|n(q)j 2 ) = q 2 G(q) obtained from atomistic 
MC simulations for the case of tension (solid lines), and com- 
pression (dashed lines), for different values of external strain. 



the range of accessible wave-vectors with respect to the 
unstrained case. In Fig. [4] we compare the correlation 
functions obtained from numerical simulations (full lines) 
to the SCSA results (dashed lines) for different values 
of the strain. To highlight the change of slope of G(q) 
due to strain, in Fig. [4] we plot q 2 G(q) tha t be comes 

First, 



III 



fiat when G(q) oc q~ as discussed in Sec. 
one notices that the G(q) calculated by atomistic sim- 
ulations deviates from those calculated in the contin- 
uum limit for wave-vectors close to the Bragg peak at 
q = |^ = 2.94 A" 1 with a = 1.42 A being the carbon- 
carbon distance in graphene. We mention also that, for 
the strained cases, the error bars of the Monte Carlo sim- 
ulations are negligible. 

Starting from the unstrained case (u = 0), we see that 
the MC (full black line) and the SCSA (dashed black 
line) results agree reasonably well and both deviate from 
the correlation function in the harmonic approximation 
(dot-dashed black line) at small wave-vectors, pointing 
out the importance of anharmonic effects at long scales in 
unstrained samples.^ For the strained cases, the SCSA 
and MC results are also comparable, what justify the 
use of SCSA when dealing with samples under tension. 
However, we must emphasize that for 0.4% strain, there 
is almost no difference between Go(q) [as obtained by Eq. 
^] and G(q) in the SCSA [full solution of Eq. @-i[l4|)], 
as discussed in Sec. 



|III| and they are exactly the same 
for the highest value of strain shown here, 1.5%. A more 
rigorous check of the validity of SCSA would require MC 
simulations for samples under even weaker strain, which 
requires several times larger samples to achieve the same 
accuracy and makes such simulations much more time 
consuming. Nevertheless, the present results already con- 
firm that rather weak strain is enough to suppress the an- 
harmonicities in stiff membranes as graphene, as it can be 



seen in Fig. [4] by the almost fiat line-shape of g 2 G(q) as 
we move from the Bragg peak towards small wave- vectors 
of the spectrum, for tensions > 0.4%. 

The study of a compressed membrane is more deli- 
cate because of the fact that its equilibrium state does 
not correspond to the flat phase any more! 23 " 26 ! There- 
fore, the standard elastic theory that we have used in 
Sec. |III| does not apply to this case. However, one can 
at least study the system by means of Monte Carlo sim- 
ulations. In Fig. [5] we compare the MC results for the 
correlation function of a tensioned membrane to that of 
a compressed membrane. There we see that G(q) for a 
compressed membrane does not show the characteristic 
gradual crossover to another power law as in the ten- 
sioned case. Instead, there exist a wave-vector for which 
the correlation function suffers an abrupt deviation from 
the harmonic behavior, presenting signatures of a possi- 
ble first order buckling phase transition. 

Finally, all our results are graphically summarized in 
Fig. |6j where we show a snapshot of the Monte Carlo 
sample for a tensioned graphene membrane [Fig. [6]ja)], 
for an unstrained membrane, Fig. [6](b) , and for a com- 
pressed graphene sheet, Fig. |6jc). In the first case, the 
equilibrium state is an almost perfectly flat membrane, 
for which the anharmonic effects have been suppressed 
due to the application of tension. The anharmonic cou- 
pling between bending and stretching modes is instead 
important for the case of an unstrained membrane, as 
the one of Fig. [6jb) , which leads to a corrugated low 
energy phase due to the existence of thermal ripples in 
the system pi and which is well described by means of 
a continuum elastic theory as the SCSA. However, this 
theory is not applicable to a compressed membrane as 
the one shown in Fig. [6^c), for which the sheet buckles 
into shapes that remove in-plane compression, in order 
to reduce its elastic energyPS 



FIG. 6. (Color online) Typical Monte Carlo configurations of a graphene sample of 8640 atoms for: a) 1.5% tension, b) 
unstrained, and c) 1.5% compression. 



V. CONCLUSIONS 



In summary, we have studied the effect of external 
strain in the correlation function of fiexural modes in the 
SCSA. In the presence of strain, three different regimes 
can be distinguished in the dispersion relation of fiexu- 
ral phonons: ojyj(q) ~ q in the long wavelength limit, 
Ufiio) ~ g 2 ~''/ 2 in the intermediate range of wave- 
vectors of the spectrum (where r\ sa 0.82 is a charac- 
teristic exponent!^ ) and finally w/j(q) ~ q 2 at shorter 
wavelengths. The results show that, for a soft mem- 
brane, rather high values of strain are needed to sup- 
press anharmonic effects, whereas for a stiff membrane as 
graphene, anharmonic effects are completely suppressed 
by less than 1% tensile strain. The correlation functions 
obtained with the SCSA compare well with those cal- 



culated from atomistic MC simulations. Taking into ac- 
count that the scattering of electrons by fiexural phonons 
has been shown to be the main limitation for the charge 
mobility in suspended graphenef^ our results point that 
the application of a small tension to the graphene layer 
would reduce the out-of-plane vibrations that lead to the 
fiexural modes, increasing the mobility of the suspended 
samples. 
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